Traction image method for irregular free surface boundaries in finite difference seismic wave simulation   总被引：16，自引：1，他引：16  

Wei Zhang  Xiaofei Chen 《Geophysical Journal International》2006,167(1):337-353

In this study, we propose a new numerical method, named as Traction Image method, to accurately and efficiently implement the traction-free boundary conditions in finite difference simulation in the presence of surface topography. In this algorithm, the computational domain is discretized by boundary-conforming grids, in which the irregular surface is transformed into a 'flat' surface in computational space. Thus, the artefact of staircase approximation to arbitrarily irregular surface can be avoided. Such boundary-conforming gridding is equivalent to a curvilinear coordinate system, in which the first-order partial differential velocity-stress equations are numerically updated by an optimized high-order non-staggered finite difference scheme, that is, DRP/opt MacCormack scheme. To satisfy the free surface boundary conditions, we extend the Stress Image method for planar surface to Traction Image method for arbitrarily irregular surface by antisymmetrically setting the values of normal traction on the grid points above the free surface. This Traction Image method can be efficiently implemented. To validate this new method, we perform numerical tests to several complex models by comparing our results with those computed by other independent accurate methods. Although some of the testing examples have extremely sloped topography, all tested results show an excellent agreement between our results and those from the reference solutions, confirming the validity of our method for modelling seismic waves in the heterogeneous media with arbitrary shape topography. Numerical tests also demonstrate the efficiency of this method. We find about 10 grid points per shortest wavelength is enough to maintain the global accuracy of the simulation. Although the current study is for 2-D P-SV problem, it can be easily extended to 3-D problem.  相似文献   

On the effect of temperature and strain-rate dependent viscosity on global mantle flow, net rotation, and plate-driving forces     

Thorsten W. Becker 《Geophysical Journal International》2006,167(2):943-957

  相似文献   

Verification of the motion fields over Southeast Asia by the Florida State University (FSU) superensemble forecast     

Mastura Mahmud  Robert S. Ross 《Singapore journal of tropical geography》2006,27(2):176-190

The Florida State University (FSU) multimodel superensemble forecast is evaluated against several other operational weather models for the Southeast Asia region. The superensemble technique has demonstrated its exceptional skills in forecasting precipitation, motion and mass fields compared to either individual global operational or ensemble mean forecasts. The motion field investigation for the season of 2001 reveals that the superensemble forecasts are closer to the observed data compared to the other global member operational models through its low systematic errors at the 850 hPa level. The FSU multimodel superensemble forecasts exhibit the lowest root mean square errors (RSMEs), the highest correlation against the best observed data and the lowest systematic errors compared to the other operational model members. These forecasts have the potential to provide better daily weather predictions over the Southeast Asia region, particularly during the early northeast monsoon that often causes heavy rainfall in the equatorial part of the Southeast Asia region.  相似文献   

Effect of lithospheric root on decompression melting in plume–lithosphere interaction models     

A. Manglik  U. R. Christensen 《Geophysical Journal International》2006,164(1):259-270

  相似文献   

The Dahuiyeh (Zarand) earthquake of 2005 February 22 in central Iran: reactivation of an intramountain reverse fault   总被引：3，自引：0，他引：3  

M. Talebian  J. Biggs  M. Bolourchi  A. Copley  A. Ghassemi  M. Ghorashi  J. Hollingsworth  J. Jackson  E. Nissen  B. Oveisi  B. Parsons  K. Priestley  A. Saiidi 《Geophysical Journal International》2006,164(1):137-148

  相似文献   

Surface wave tomography: finite-frequency effects lost in the null space   总被引：2，自引：1，他引：2  

Jeannot Trampert  Jesper Spetzler 《Geophysical Journal International》2006,164(2):394-400

  相似文献   

Stress Green's functions for a constant slip rate on a triangular fault     

Taku Tada 《Geophysical Journal International》2006,164(3):653-669

  相似文献   

A traveltime reciprocity discrepancy in the Podvin & Lecomte time3d finite difference algorithm     

Ari Tryggvason  Björn Bergman 《Geophysical Journal International》2006,165(2):432-435

We have found that the extensively used finite difference scheme time3d produces time fields which are 'asymmetric' in the sense that traveltimes computed to the right of the source are faster than traveltimes computed to the left. All finite difference schemes will, as they are approximations to the wave equation, to some extent fail to obey reciprocity perfectly. We show, however, that the errors in time3d may be significant—and unnecessarily large. An asymmetry in the point source initialization has been identified, and after correction time3d produces time fields with an improved reciprocity.  相似文献   

An integral equation approach for 2.5-dimensional forward and inverse electromagnetic scattering     

Aria Abubakar  Peter M van den Berg  Tarek M Habashy 《Geophysical Journal International》2006,165(3):744-762

  相似文献   

A simple hypoplastic model for normally consolidated clay     

Wen-Xiong Huang  Wei Wu  De-An Sun  Scott Sloan 《Acta Geotechnica》2006,1(1):15-27

The paper presents a simple constitutive model for normally consolidated clay. A mathematical formulation, using a single tensor-valued function to define the incrementally nonlinear stress–strain relation, is proposed based on the basic concept of hypoplasticity. The structure of the tensor-valued function is determined in the light of the response envelope. Particular attention is paid towards incorporating the critical state and to the capability for capturing undrained behaviour of clayey soils. With five material parameters that can be determined easily from isotropic consolidation and triaxial compression tests, the model is shown to provide good predictions for the response of normally consolidated clay along various stress paths, including drained true triaxial tests and undrained shear tests.  相似文献